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c - Specific  heat  of  solid. 

B - Energy  absorbed  at  ignition. 

q - Average  heat  flux  defined  by  Equation  (5). 

q - Heat  flux. 

r - Regression  rate  of  solid. 

t.  - Ignition  delay, 

ig 

t - Tine. 

Tg  - Surface  temperature  of  solid. 

T - Initial  tenperature. 

o 

V/S  - Volume-to-surface  ratio  for  a grain. 

a - Thermal  dif fusivity  of  solid. 

X - Thermal  conductivity  of  solid. 

p - Solid  propellant  density. 

AT.  - Bulk  temperature  change. 
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I.  INTRODUCTION 


It  is  no  secret  that  it  is  difficult  to  calculate  the  ignition 
sequence  in  a bed  of  solid  propellant  grains.  Computer  codes  have 
attempted  to  describe  it  by  solving  one-dimensional  equations  for  the 
two  phase  flow.  Even  if  one  concedes  correctly  posed  differential 
equations,  there  remains  the  question  of  posing  the  auxiliary  equations 
or  data  for  ignition  criteria,  heat  transfer  to  the  grains,  transient 
regression  rates,  igniter  gas  flow,  and  the  like. 

Application  of  these  computer  codes  to  solid  propellsnt  guns  has 
principally  focussed  on  the  description  of  pressure  wave  development 
in  the  gun  chamber.  More  than  any  other,  it  has  become  the  standard 
by  which  to  judge  the  merit  of  a code.  And  well  it  should  because 
pressure  difference  is  one  of  the  few  presently  measurable  phenomena  to 
compare  with  prediction.  Parametric  studies  have  established  that  re- 
gression rates , bore  resistance , and  drag , can  be  varied  enough  to  force 
agreement  between  predicted  and  measured  pressures.  For  the  time  being, 
this  agreement  is  taken  as  a sufficient  basis  to  Justify  interior  ball- 
istics studies  of  gun  behavior  as  the  propelling  charge  is  modified^ 

There  is  considerable  evidence  that  the  ignition  behavior  strongly 
Influences  the  gun's  susceptibility  to  catastrophic  overpressure. 
Describing  the  ignition  requires  first  describing  how  the  propellant 
itself  makes  the  transition  from  heating  to  combustion.  In  the  gun, 
the  ignition  front  movement  is  the  progression  of  the  transition  event 
through  the  bed.  The  one  measure  of  that  movement  is  the  flamespreading 
rate  which  is  only  beginning  to  be  measured.  Such  measurements  will, 
however,  provide  another  standard  by  which  to  judge  a code.  It  is 
necessary  now  to  examine  the  details  of  the  ignition  and  its  computation 
to  improve,  if  feasible,  the  postulated  laws. 

This  report  will  address  the  ignition  criterion  for  the  propellant 
and  the  heating  of  the  solid  to  the  ignition  regime.  It  will  examine 
the  existing  experimental  ignition  data  and  computational  technique 
of  the  NOVA  computer  code  with  a view  to  improving  the  postulated  igni- 
tion sequence. 

II.  IGNITION  CRITERION 

It  would  be  most  convenient  if  the  solution  of  the  transient  heat 
conduction  in  the  solid  led  to  an  unambiguous  condition  that  marked  a 
sharp  transition  from  endothermic  heating  to  self-sustained  combustion. 
Unhappily  such  a condition  cannot  be  found  by  reference  to  experimental 
data  on  Ignition  of  nitrocellulose  base  propellants.  Easily  found  in 
the  conduction  calculation  are  surface  temperature,  total  energy  ab- 
sorbed, distribution  of  temperature  and  energy;  of  these,  none  has  been 
uniquely  correlated  with  ignition.  Worse,  the  propellant  sees  a time- 
dependent,  convective  energy  source  but  all  the  experiments  used  a 
constant  radiant  source.  It  is  nonetheless  necessary  to  somehow  apply 
the  data  to  obtain  a usable  ignition  criterion. 


total  absorbed. energy  as  a criterion,  usually  without  detailed  Justifi- 
cation. Gough1  referenced  the  theory  of  Andersen2  but  used  surface 
temperature  for  simplicity  because  "the  availability  of  supporting  data 
favors  (it)".  It  Is  rather  the  lack  of  available  data  on. the  parameters 
of  Andersen's  theory  that  precludee  Its  use.  Koo  and  Kuo  postulated  a 
sequence  In  which  ablation  started  at  some  threshold  temperature.  The 
transition  to  full  combustion  at  the  ignition  temperature  allowed  a 
convective  contribution  linearly  dependent  on  how  near  the  surface  tem- 
perature Is  to  the  full  Ignition  temperature.  This  is,  in  effect,  a 
surface  temperature  criterion  because  no  other  feature  is  considered. 
References  for  the  specific  temperatures  were  not  reported  although  they 
are  within  the  generally  accepted  limits  for  ablation  and  combustion  of 
nitrocellulose  base  propellants. 


Krier  selected  a criterion  consistent  with  his  continuum  mechanics 
approach  by  specifying  a bulk  temperature  rise  in  the  solid  phase. 
Although  this  is,  on  its  face,  a total  energy  absorbed  criterion,  it 
cannot  be  directly  compared  with  ignition  experimental  data.  Ignition 
energy  is  reported  on  a unit  surface  area  basis;  the  bulk  temperature 
criterion  is  on  a unit  mass  basis.  To  compare  the  two,  consider  Krier' s 
gun  calculation.  For  a bulk  temperature  rise  ATg  the  energy  absorbed 
per  unit  surface  area  is 


p - density 
c ” specific  heat 
(V/S)  ■ volume  to  surface  ratio 


1.  P.S.  Gough , "Numerical  Analysis  of  a Two-Phase  Plow  with  Explicit 
Internal  Boundaries",  IHCR  77-5 , Naval  Ordnance  Station,  Indian 
Head,  MD,  April  1977. 

2.  W.H.  Andersen,  "Model  of  Transient  Ignition  to  Self  Sustained 
Burning ",  Comb.  Soi.  and  Tech.,  5,  p.  75-81  (1972). 

3.  J.H.  Koo  and  K.K.  Kuo , "Transient  Combustion  in  Granular  Propellant 
Beds , Part  I.  Theoretical  Modeling  and  Numerical  Solution  of 
Transient  Combustion  Processes  in  Mobile  Granular  Propellant  Beds", 

BRL  Contract  Report  No.  346,  August  1977.  (AD  #A044998) 

4.  H.  Krier,  "Predictions  of  Flamespreading  and  Pressure  Nave  Propagation 
in  Propellant  Beds",  University  of  Illinois  Technical  Report  AAE- 
75-6,  July  1975. 
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It  should  be  noted  here  that  for  a given  Ignition  energy,  presumably 
obtainable  from  experimental  data,  the  bulk  temperature  assigned  depends 
on  the  volume  to  surface  ratio  of  the  grains.  In  the  problem  demonstrate 
by  Krier,  a 76mm  gun,  the  ratio  for  the  seven  perforated  grains  was 
0.014  cmJ/cm2.  If  the  gun  contained  ball  propellant  the  ratio  mould 
have  been  about  0.18  for  balls  having  the  same  mass  as  the  perforated 
cylinders.  With  such  a spread  in  the  geometry  term,  there  would  have 
to  be  a corresponding  spread  of  an  order  of  magnitude  for  the  bulk 


temperature  rise  to  keep  the  ignition  energy  fixed.  In  the  76asn  gun 
simulation,  Krier  used  a nominal  bulk  temperature  rise  of  30  R (17K) 
which  presumes  an  ignition  energy  of  0.46  3/ car-. 

When  equivalent  time  to  ignition  was  the  criterion  for  evaluating 
the  bulk  temperature  concept,  a different  conclusion  was  reached  earlier 
For  a constant  heating  rate  of  750  W/cm2  a bulk  temperature  rise  of  4K 
gives  the  same  ignition  time , as  a surface  temperature  rise  of  155K  for 
a multi-perforated  grain.  This  should  not  be  considered  an  indictment 
of  the  bulk  temperature  concept  any  more  than  the  surface  temperature 
concept.  Both  have  absorbed  the  same  energy  at  ignition,  0.17  J/cm  , 
a value  low  by  reference  to  any  experimental  data. 

Another  popular  criterion  is  that  the  ignition  energy  must  approxi- 
mate the  energy  content  of  the  solid  for  steady-state  combustion  at  the 
ignition  pressure.  For  composite  propellants,  Derr  and  Fleming”  found 
a correlation  of  ignition  delay  with  steady  state  burning  rate  for  any 
value  of  the  constant  stimulus.  If  one  assumes  (1)  burning  rate  to  be 
at  function  of  surface  temperature  only  and  (2)  no  subsurface  reactions, 
the  energy  in  the  solid  is 


Tg  - surface  temperature 
a - thermal  diffuslvity 

,.T  ' 

r " regression  rate 


C.W.  Nelson,  "Comparison  of  Predictions  for  Three  Two-Phase  Flow 
Codes ",  BRL  Memorandum  Report  No.  2729,  February  1977  (See  also 
13th  JANNAF  Combustion  Meeting,  September  1976).  (AD  #A0S7348) 

R.L.  Derr  and  R.W.  Fleming,  "A  Correlation  of  Solid  Propellant  Arc 
Image  Ignition  Data",  10th  JANNAF  Combustion  Meeting,  Newport,  RI, 
CPIA  Publication  243,  August  1973. 


At  a given  burning  rata,  tha  energy  should  than  ba  independent  of  tha 
flux  level;  for  tha  0.5  ca/sec  burning  rata,  tha  anargy  la  about  5.4 
I /car.  Proa  tha  axparlaantal  data 

q (W/ca2)  420  210  125  63 
E (J/ca2)  2.7  3.8  3.8  4.5 

Total  energy  (E)  is  neither  constant  nor  simply  dependent  on  heat  flux 
(q).  The  steady  burning  energy  is  greater  than  the  ignition  energy  by 
at  least  30Z. 


For  nitrocellulose  propellants  with  the  values  used  in  earlier 
studies 7,  the  steady  staite  energy  in  the  solid  at  a regression  rate  2 
0.75  ca/sec,  is  only  about  1.3  J/cm2.  The  DeLuca  data*  reported  1.9  J/cin 
energy  at  210  tf/cm2  flux  and  .84  J/cm2/sec  at  1050  W/cm2  (extrapolated). 
These  values  are  independent  of  pressure  and  therefore  burning  rate. 

It  is  then  not  obvious  how  steady  burning  energy  content  and 
ignition  are  directly  related. 


NOVA  code  studies  of  the  155mm  howitzer  with  M30  propellant  have 
to  date  assumed  a surface  temperature  ignition  criterion  of  450K,  a 
value  obtained  by  using  the  DeLuca**  data  for  the  arc  image  ignition  of 
the  highly  carbonized  double  base  propellant  with  a slope  of  -2.0, 
that  is. 


tig«q 


-2.0 


Jr;,: 


where 


ignition  delay  time 
heat  flux. 


Data  on  constant  radiant  flux  stimulation  of  nitrocellulose  pro- 
pellants are  available  from  the  experiments  of  DeLuca  et  al&  for  both 
arc  image  and  laser  sources.  They  found  that  some  propellants  exhibit 
a slope  of  -2.0  which  means  Ignition  energy  and  ignition  surface  tem- 
perature are  independent  of  heating  rate.  In  most  cases  the  relationship 
was  t a q~l*6  which  means  some  mechanisms  other  than  simple  surface 
heat  absorption,  presumably  gas  phase  chemistry  or  in-depth  absorption 
and  heating.  Because  some  of  the  propellant  formulations  absorb  arc 


7.  L.H.  Caveny , M.  Swmerfield , and  C.W.  Nelson,  ", Ignition  Transients 
and  Pressurization  in  Closed  Chambers ",  BRL  Memorandum  Report  No. 
2558,  November  1975.  (AD  #A017747) 

8.  L.  DeLuoa,  L.H.  Caveny,  T.J.  Ohlemiller,  and  M.  Surmerfield, 

" Radiative  Ignition  of  Double  Base  Propellants ",  Part  I.,  14_  (7), 
p.  940-946 , Part  II,  14  (8),  pp.  1111-1117  AIAAJ  (1976). 
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radiation  In— depth  and  because  the  propellant  stimulus  in  the  gun  chamber 
is  probably  convective,  it  is  reasonable  to  consider  only  the  laser 
heating  data  which  is  practically  all  absorbed  at  the  surface  as  con- 
vective flux  would  be. 

Corroboration  of  these  Deluca  data  would  be  expected  from  arc  image 
experiments  of  Hayes  and  Lenchltz^  who  studied  ignition  of  standard 
Army  propellants.  But  they  found  slopes  of  -1.0  to  -1.3. 

Table  I and  Figure  1 compare  the  two  sets  of  experimental  data  for 
uncatalyzed  nitrocellulose  propellants.  The  heat  flux  to  the 
propellant  surfaces  in  the  granular  bed  has  been  computed  (by  NOVA  code) 
to  grow  roughly  linearly  in  time  to  about  1250  W/cm^  at  ignition.  A 
column  has  been  added  to  predict  the  ignition  time  at  a constant  flux  of 
840  W/cm2  (an  arbitrary  but  reasonable  "average”  value)  assuming  that 
the  time  vs  flux  data  can  be  extrapolated  with  no  change  in  slope,  from 
a reference  value  of  210  W/cm^,  i.e., 

tlg(S4°)  - tlg  (210)  («)" 
where  n is  the  reported  slope. 


; 


The  specific  criterion  is  needed  in  the  155mm  howitzer  application 
for  the  triple  base  M30  propellant.  Bayes  and  Lenchitz  Ignited  both  a 
usual  film  sample  and  a cylindrical  grain  in  air  with  both  a laser  and 
an  arc.  The  laser  fluxes  were  too  low  to  be  applicable  and  the  results 
too  different  from  arc  image  results  to  be  credible. 

For  M26  (double  base)  and  Ml  (single  base)  propellants,  one  should 
expect  comparability  with  DeLuca's  findings  for  Plastisol  (double  base) 
and  M9  (single  base)  propellants  respectively.  Unfortunately,  the 
comparison  is  weak.  The  slope  difference  for  the  double  base  is  —1.5 
vs  -1.0,  and  for  the  single  base  -2.0  vs  -1.0. 

Russian  radiant  flux  data  reports  ignition  of  nitrocellulose  base 
propellant  at  low  flux.  Over  the  range  5-40  W/cm^  the  slope  was  -1.6, 
consistent  with  DeLuca,  but  of  only  limited  supporting  value. 

Given  that  this  data  must  somehow  yield  an  estimate  of  ignition 
time  vs  flux  level,  it  seems  that  the  most  likely  candidate  is  the 
data  of  DeLuca  for  the  laser  heating  in  air  with  a slope  of  -1.6  and 
a predicted  delay  of  1.0  ms  at  840  W/cm^  flux. 


9.  E.  Hayee  and  C.  Lenchitz,  "The  Thermochemistry  of  Ignition",  BRL 
Report  No.  183 3,  ed.  E.  Freedman,  K.J.  White,  I.W.  May,  1976. 

(AD  #A0 16902) 
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TABLE  I.  NC  BASE  PROPELLANTS  IGNITION 


For  convective  hotting  to  ignition,  Keller  et  al10  used  a shock 
tunnel  with  a constant  gaa  velocity  inlet  gas  temperature  to  IgnitOj 
composite  propellant  samples  using  heat  fluxes  from  125  to  630  tf/cm  . 

The  slope  was  -2.0  compared  to  -1.8  in  the  radiation  fumaca.  Grigorev 
at  al11  ignited  pyroxylin  (nltrocelluloee)  cylinders  by  alow  convective 
heating  (10  sec)  [convective  heating  coefficient,  h * .01  W/cm?/*K  ,, 
compered  to  a calculated  h « 0.4  in  the  gun  Chamber] . Kosenband  at  alx 
ignited  pyroxylin  slabs  with  combustion  product  gases  with  moan  heat 
fluxes  8-50  W/cm2  although  in  each  test  the  flux  was  monotoolcelly  in- 
creasing. Plotting  ignition  time  vs  mean  heat  flux  yields  a slope  of 
-3  with  pyroxylin  products  and  -1.5  to  -2.0  with  composite  propellant 
products.  A scatter  of  60Z  in  ignition  timed  was  reported.  Rorenband 
at  all 3 reported  that  in  an  incraaslng  heat  flux  the  effect  of  [clastic 
and  thermophysical  parameters  on  ignition  time  is  much  lea*  than  in  a 
decreasing  or  constant  heat  flux.  For  rapidly  increasing  beat  flux, 
the  effect  "may  be  unimportant  and  ignition  lag  can  be  determined  solely 
by  the  dynamics  of  the  external  heat  flow".  Marshanov  and  Averson}-*  in 
a review  paper  observed  that  thermal  theory  and  experiment  agree  to 
within  about  10Z.  They  also  observed  that  in  increasing  heat  flow, 
ignition  times  are  nearly  the  same. 

Hot  wire  ignition  yields  no  more  credible  data.  Zarlco  and  Khlevnoi^ 
ignited  balllstite  (double  base)  propellant  with  an  imbedded  nickel  wire 
using  mean  heat  fluxes  ranging  from  20-125  W/cm2.  Although  experimental 
results  "agreed"  with  their  theoretical  calculations,  the  slope  of  the 
ignition  time  ys  mean  heat  flux  was  near  -4.0  for  the  experiment  and 
-3.0  for  the  theory.  Heat  flux  history  was  net  reported.  Accuracy  of 
total  ignition  energy  was  estimated  to  be  not  less  than  20Z. 


10.  J.A.  Keller,  A.D.  Baer , and  N.W,  Ryan , "Ignition  of  AP  Composite 
Propellants  by  Convective  Heating",  AIAA  J,  £,  (8),  p.  1358-1365 
(1966). 

11.  I.M.  Grigorev,  V.I.  Lisitekii,  and  A.G.  Merzhanov , " Ignition  of 
Condensed  Substance  in  Hot  Gas",  Fizika  Goreniya  i Vzryva,  3_  (4), 
p.  512-526  (1967). 

12.  V.I.  Rozenband,  V.V.  Borzykin , and  A.G.  Merzhanov , "Ignition  of 
Condensed  Substances  by  Convective  Heat  Fluxes  of  Medium  Intensity 
Under  Dynamic  Conditions ",  Fizika  Goreniya  i Vzryva,  £ (2),  p. 
171-175  (1968). 

! 

13.  V.I.  Rozenband,  A.E.  Averson,  V.N.  Borzykin,  and  A.G.  Merzhanov, 
"Some  Characteristics  of  Dynamic  Ignition  Regimes ",  Fizika  Goreniya 
i Vzryva,  £ (4),  p.  494-500  (1968). 

14.  A.G.  Merzhanov  and  A.E.  Averson,  "The  Present  State  of  Thermal 
Ignition  Theory ",  Combustion  and  Flame,  16,  p.  89-124  (1971). 

15.  V.E.  Zarko  and  S.S.  Khlevnoi , "Hot  Wire  Ignition  of  Ballistite 
Propellants",  Fizika  Goreniya  i Vzryva,  £ (2),  p.  158-170  (1968). 


Given  the  data  acattar  and  the  difficulty  of  interpreting  the 
slopes  of  >3  with  little  or  no  interned late  information,  there  feena 
little  purpose  in  trying  to  apply  it  to  the  gun  ignition  aeqaance . A 
more  useful  courae  a sens  to  be  to  adapt  the  radiant  flux  experimental 
data  to  the  unsteady  convective  heating  la  the  gun* 

The  Princeton  laser-in-alr  data  yields  a formula 

tig  " 3*98  • <3) 

where  q is  the  Met  flux  in  cal/cn2/sec;  tlg  is  in  sac.,, 

A typical  NOVA  calculation  shows  the  heat  flux  to  a grain  to  be  roughly 

-Oi  ' r j-'i  t *1  XV  H 'Mahil  .1  r 2©  ■ Z 1 Yi&u&ltfiit'l  gf/s> 

q - 13.4  + 4.6  x 105(t-to)  , 

where  the  heating  starts  at  tQ. 

The  mean  heat  flux  is  then 


q = to 

o to 


t 

/ qdt 


One  route  to  computing  the  ignition  delay  would  be  tha  expression 

tig  - s.psc^r1*6  , (5) 

•*  ' f.  • ■ : • ll  aJiiK.K  i 

although  this  may  not  work  well  for  a rapidly  changing  flux  such  as 
tlg  - *.98  [13.4  + 4,6  *- tig]"1,6  . <*) 


Solve  to  find 


t . - 0.8  ms, 

lg 


■ n%iioi 

jii <13  t i 


whereas  the  surface  temperature  criteria  predicts  about  0.3 -ns. 

- . ‘j  . ■ • .7  C 1 t 'j:i  : ■ - / • : $.:•  fBtfi 

An  alternative  i^  to  use  the  average  value  of  (q”sl,-> 


(q‘1,6> 


avg  t-t 


/*’  q-1*6^ 


which  for 


q - a + b (t-to) 


Whan  this  expression  Is  substituted  into  Bq.  (4),  Which  is  than  solved 
for  Ignition  tine 


This  is  s long  tine  compared  to  0.3  ns. 

The  heat  flux  estimates  were  token  froa  calculations  Which  only 
went  to  the  0.3  ns  ignition.  The  0.8  ns  estimate  depends  on  extra- 
polation of  the  heat  flux.  One  eheck  Is  to  calculate  the  flanespresding 
using  the  continuously  calculated  mean  heat  flux  to  decide  if  ignition 
criterion  has  been  satisfied.  Eq.  (4)  becoaea 


If  t > (t.  + t ) , ignition  has  occurred 


A test  case  for  these  ignition  studies  can  be  constructed  easily 
by  modifying  the  input  data  for  a 155ns  howitzer  calculation.  The  only 
changes  are  shown  in  Table  II. 

TABLE  II.  INPUT  DATA  CHANGES  FOR  TEST  CASE  (NOVA  CODE) 


Test  Case 


Movable 
82  ca 
53  cm 
5 ca 


Projectile 
Chamber  Length 
Charge  Length 
Igniter  Length 

These  changes  are  not  enough  to  affect  the  flaaespreadlng  in  aost  of 
the  bed,  but  they  do  shorten  the  calculation. 

Using  the  ignition  delay  criterion  (-1.6  slope)  reduces  t]*e  calcul- 
ated flaaespreadlng  rate  markedly.  Whereas  the  rate  was  about  6.27  aa/i 
for  the  fixed  ignition  temperature  (450  K).  it  is  only  half  that  value 
with  the  average  flux  method.  Froa  a rate  of  0.1  ns/ pa  just  outside 
the  igniter  input  region  it  reaches  0.15  ns/ue  in  a distance  qf  17  ca; 
thereafter  it  is  approximately  constant.  The  predicted  propagation 
rate  of  the  first  absorption  of  heat  by  the  grain  decelerates  froa 
Q.6  ns/vs  near  the  Igniter  to  0.35  ns/ps  downstream.  (The  sound  speed 
in  the  undisturbed  air  in  the  bed  is  0.34  aa/ps.)  Ignition  delay  time 
Increases  from  0.4  as  at  the  igniter  region  to  1.9  as  well  downstreaa. 


Ml  these  treads  can  be  seen  in  Figure  2 where  they  ere  compared  with 
the  ignition  speed  predicted  with  the  fixed  surface  temperature 
criteria  of  810  and  1500  R. 

This  ignition  criterion  end  transient  regression  cannot  be  simul- 
taneously applied.  The  transient  regression  assumes  that  the  temperature 
profile  in  the  solid  at  ignition  is  the  initial  condition  for  the 
combustion.  This  ignition  criterion  Ignores  the  temperature  profile 
and  looks  only  at  the  rate  of  heat  absorption.  The  gasification  period 
of  the  flame  development  is  not  a straightforward  conduction  problem 
And  the  temperature  profile  cannot  be  so  computed.  Until  a nodal  for 
the  flame  development  can  give  initial  conditions  for  the  combustion, 
a quasi-steady  pressure  dependence  is  the  only  credible  formula. 

111.  INTEGRAL  vs  FINITE  DIFFERENCE  SOLUTION 

Heating  of  the  solid  to  ignition  requires  a solution  to  the 
transient  heat  conduction  equation.  Because  the  boundary  condition  of 
heat  flux  cannot  be  specified,  a numerical  solution  is  dictated.  A 
choice  must  be  made  between  "exact"  finite  difference  and  approximate 
techniques.  Both  have  been  reasonably  well  studied;  it  remains  only 
to  judge  their  relative  errors  for  the  problem  at  hand. 

It  is  widely  recognized  that  the  FDE  solutions  can  com  almost 
as  close  to  the  exact  solution  as  the  user  is  willing  to  shrink  his 
discretization.  Kooker  and  Nelson1^  found  that  for  an  intense  heating 
problem  the  invariant  imbedding  technique  estimated  surface  temperature 
within  1%  of  the  exact  solution. 

A well  known  approximate  method  for  transient  conduction  is  the 
Integral  method  of  Goodman  which  has  a surface  temperature  error  of  9% 
for  a parabolic  and  2%  for  a cubic  profile  when  the  Incident  heat  flux 
is  constantly.  Because  the  heat  flux  in  the  propellant  bed  is  computed 
to  be  approximately  a linear  ramp  in  time,  it  would  be  worthwhile  to 
estimate  the  error  for  these  profiles. 
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The  form  of  the  exact  solution  is  given  by  Carslaw  and  Jaeger  as: 


16.  D.E.  Kooker  and  C.W.  Nelson , " Numerical  Solution  of  Solid  Propellant 
Transient  Combustion ",  ASMS  Paper  77-HT-17 , August  1977. 

17.  T.R.  Goodman , "The  Heating  of  Slabs  with  Arbitrary  Heat  Inputs ", 

J.  Aerospace  Sciences , 187-188,  March  1969. 

18.  H.S.  Carslaw  and  J.C.  Jaeger,  "Conduction  of  Heat  in  Solids", 

Oxford  Press,  2nd  ed.,  pp.  78,  1969. 
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Flame  Spreading  Comparisons  of  Ignition  Criteria 


For  the  cubic  profile  where 


/ q(t)dt 


ho 11c  profile  will  predict  e higher  surface 
4 or  about  6Z  of  the  temperature  rise. 


Comparison  with  the  exact  solution  can  be  obtained  by  using  the 
typical  gun  heating  rate 

q - 56  + 1.9  x 106t  (W/cm2)  (13 


which  estimates  an  error  of  15Z  for  the  parabolic  profile  and  8Z  for 
the  cubic  profile.  For  the  more  useful  estimate  of  error  in  time  to 
reach  a specified  surface  temperature  the  eetlmates  are  10Z  for  the 
parabolic  and  6Z  for  the  cubic. 

The  FDE  Invariant  imbedding  technique  produces  en  estimated  error 
of  leee  than  1Z. 

Since  there  le  little  difference  in  computational  effort  between 
the  parabolic  and  cubic  approximations  it  is  not  surprising  that  the 
cubic  is  already  used  by  Gough  and  Kuo.  And  if  surface  temperature 
accuracy  were  the  only  criterion,  it  would  be  rather  straight  forward 
to  make  the  trade-off  between  accuracy  and  efficiency.  If  quasi-steady 
combustion  is  assumed  after  ignition,  eurface  temperature  is  probably 
an  adequate  teet  because  the  Ignition  criterion  will  probably  be 
dominated  by  the  surface  behavior. 

! + J?  ’’  P 

If.  however,  transient  combustion  is  to  be  calculated  with  a 
thermal  theory  model,  it  is  also  necessary  to  have  an  accurate  temperature 
distance  profile  in  the  solid.  For  this  reason,  the  only  code  coupling 
tranelent  combustion  to  the  flow  uses  an  FDE  technique  throughout.  For 
a constant  heating  problem,  the  error  in  the  temperature  vs  distance 
calculation  can  be  found  by  direct  calculation  of  both  the  exact  and 
the  integral  solutions.  The  error  of  2.3Z  at  the  surface  represents 
a temperature  error  of  8.3  K for  a problem  of  beating  at  840  W/cm2  for 
0.8  me.  The  temperature  error  decreases  with  distance  to  a distance 
of  1.5VotT  at  which  point  the  two  solutions  cross  over  such  that  the 
approximate  solution  is  less  than  the  exact  solution  and  the  error 
begins  to  grow  although  the  values  of  both  are  small. 

IV.  NUMERICAL  PROBLEMS 

NOVA  economizes  on  storage  and  running  time  by  using  marker 
particles  (fewer  than  grid  points)  to  compute  the  ignition  and  transient 
combustion.  Matching  the  solutions  there  to  the  flow  field  at  the  grid 
mesh  points  requires  interpolation  in  both  directions.  In  the  limit  as 
the  marker  particle  spacing  approaches  the  mesh  spacing.  Interpolation 
errors  shrink  toward  zero. 

v -jT  *•  l j ■>  T >■ 

The  original  scheme  computed  grain  temperatures  in  the  flow  meeh 
by  linear  interpolation  between  the  markers.  Surface  temperature 
computation  of  the  marker  stopped  at  ignition.  The  result  was  a leap 
frog  flamespreading  as  shown  in  Figure  3.  The  problem  was  the  tempera- 
ture Interpolation.  An  artifice  was  then  Incorporated  which  continued 
to  heat  the  ignited  marker  as  an  inert  body  without  taking  the  heat 
from  the  gas.  It  was  hoped  thereby  to  smooth  the  interpolation  without 
affecting  the  energy  balance  in  the  gas.  Things  improved  considerably 
as  ehown  by  the  curve  in  Figure  3 marked  3.0.  By  comparison  the  grid 
spacing  was  0.6,  leaving  five  (5)  grid  points  between  marker  particles. 
Even  this  flamsspr ending  is  ragged;  the  local  flamespreading  rate  varies 
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DISTANCE  (cm) 

Effect  of  Marker  Spacing 


over  a factor  of  two.  With  tha  aarkar  spacing  equal  to  the  grid  spacing 
over  a portion  of  the  bad,  tha  flaasspreadlng  rata  was  constant  (vari- 
ation laas  than  22) . Note  that  at  tha  and  of  tha  closely  packed  aarkar 
particles  the  rata  changes  draaatlcally  aa  tha  Interpolation  Interval 
jimp a to  9.5  between  there  and  the  forward  end  of  the  bad.  A clear 
route  la  thus  marked  to  a saooth  computation;  whether  the  convenient 
artifice  of  superheating  fatally  distorts  tha  physics  raaalna  to  be 
determined  by  the  effects  on  the  chamber  dynamics.  Apparently  tha 
overall  rata  Is  only  slightly  affected. 

; ' -v  / d' 

By  good  fortune  this  error  efficiency  trade-off  Is  unnecessary 
because,  as  Is  shown  in  the  section  on  FDE  vs  integral  solution,  errors 
are  small  enough  to  avoid  a full  numerical  integration  up  to  Ignition. 
The  Integral  approximation  is  adequate  and  can  be  economically  used  at 
every  grid  point. 


V.  EFFECT  OF  IGNITER 

Flamespreading  calculations  usually  show  a dependence  of  the  rate 
on  the  assumed  strength  of  the  igniter  gas  flow.  In  the  limit  of  a 
long  bed,  the  effect  should  die  out  far  from  the  gas  source  if  the 
rate  becomfs  invariant  In  time.  To  test  the  situation,  a problem  was 
created  to  assume  an  igniter  of  variable  strength  at  the  end.  Specifi- 
cally, the  igniter  occupied  the  first  3.^  cm  of  a bed  S3  cm  long. 

Table  III  shows  the  results. 
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TABLE  III.  EFFECT  OF  IGNITER  FLOW  RATE 


Assumed  Primer 

Flags 

espread  Rate  @ 

Flow  Rate  (kg/ cm/ sec) 

15  cm(WHs) 

30  cm 

40  cm 

50  cm 

A? 

.45 

.04 

.10 

.10 

.07 

.89 

.08 

.11 

.11 

.11 

1.79 

.09 

.14 

.14 

.10 

3;58 

.13 

.17 

.19 

.14 

» ; j / //  i 

Igniter  flow  rate  can  be  increased  by  changing  its  assumed  length 
while  keeping  the  flow  rate  per  unit  length  the  same.  The  effect  of 
such  a variation  is  shown  In  Table  IV  for  the  nominal  case  of  a delay 
time  criterion  and  0.9  kg/ cm/ sec  flow  rate. 

TABLE  IV.  EFFECT  OF  IGNITER  LENGTH 


Igniter  Length  (cm) 

Flame  Spread  Rate 
(mm/ys) 


3.1 

.06+. 11 


lai 

.08+. 12 


15»1. 
.10+. 15 
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®#eWM  I*  «•  usually  not  long  compared  to  Ignitor 

length,  any  calculatlou  of  fl  smeaprasding  rata  will  be  sens  It  Ira  to 
igniter  assamptloma.  Such  unpleasant  news  mane  that  assumptions 
about  the  igniter  nest  either  be  carefully  formulated  from  experimental 
data  or  the  flow  rate  will  become  just  smother  arbitrary  parameter. 

That  la  not  to  any,  however,  that  the  gun  calculation  will  be  otherwise 
useless;  soma  Important  ballistics  variables  seam  Insensitive  to 
flamespreadlng  rats.  How  Insensitive  Is  not  Immediately  assessable. 

VI.  COMPARISON  WITH  MEASUREMENTS 
19 

East  at  al  have  made  the  few  measurements  of  the  flamespreadlng 
rates  In  naval  gun  propellant  beds.  Their  data  are  shown  in  Table  ▼ 
which  shows  the  major  variables  expected  to  affect  the  flamespreadlns 
rate.  ^ 

ZABLE  V.  MEASURED  FLAMESPREADIMG  RATES 


Nominal 

Grain 

Burning 

Flamespread 

Chamber 

Length 

Diameter 

Rate  7 MPa 

Rate 

Diameter 

Propellant 

76  mm 

M6 

17 

7.4 

13 

.36-. 48 

76  mm 

M6+2 

17 

7.4 

12 

.35-. 39 

127  mm 

M26 

30 

14 

11 

.44 

127  mm 

NACO 

13 

6.9 

4 

.30 

It  may  well  be  noted  that  there  are  experimental  uncertainties  in 
igniter  gas  flow,  in  flame  front  location  at  only  a few  discrete 
locations,  and  in  distinguishing  the  ignition  front  from  the  luminous 
gas  penetration  ahead  of  the  front.  Nevertheless,  it  Is  clear  that 
the  observed  rates  are  much  higher  (0.4  vs  0.2)  than  those  calculated 
under  the  ignition  delay  criterion.  The  only  calculations  which  estimate 
these  observed  rates  are  those  under  a surface  temperature  criterion  of 
about  150#C. 

For  the  155mm  howitzer  propellant  bed  which  formed  the  basis  for 
these  calculations,  the  propellant  geometry  and  burning  rates  are 
the  same  as  the  127mm  gun.  One  would  expect  the  flamespreadlng  rate 
to  vary  only  in  relation  to  the  degree  by  which  the  Igniter  gas  flow 
differs. 


J.L.  East  and  D.R.  MoClure,  "F  > Studies  of  Ignition  and 

Combustion  in  Basal  Ouns'\  Pr^m^ngs  of  the  12th  JANNAF  Combustion 
Meeting , Beuport , RIt  August  1975  (See  also  11th  JANNAF  Combustion 
Meeting) . 
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TABLE  VII.  EFFECT  ON  BALLISTICS 


Commodity 

Experimental 

Base  Case 

Revised 

Ignition 

Doubled 

Revised 

Ignition 

Flame  Spread  Rate 

unknown 

0.35-0.44 

0.2-0. 3 

.17-. 21 

(mm/ vs) 

Peak  Pressure  (MPa) 

235 

293 

268 

265 

Pressure  Differences 

-40/35 

-73/48 

-75/51 

-67/53 

(MPa)* 

Muzsle  Velocity  (m/s) 

790 

768 

NC 

* - First  number  is  first  negative  difference;  second  is  second  positive. 

An  immediate  conclusion  drawn  from  Table  VII  is  that  the  ignition 
criterion  does  not  appreciably  affect  pressure  difference  amplitudes. 

The  results  also  show  that  the  reduction  of  50Z  in  flame  speed  reduces 
peak  pressure  by  only  10Z. 

' 

It  should  not  be  considered  significant  in  this  study  that  all 
the  simulations  over-predict  the  actual  ballistics.  This  condition 
could  be  remedied  by  adjusting  the  mostly  guessed-at  bore  resistance 
function.  The  trends  should  remain  the  same. 

VIII.  CONCLUSIONS 

(1)  Laser-in-air  radiant  heating  data  yield  the  most  credible 
basis  for  an  ignition  criterion. 

(2)  An  integral  method  based  on  a cubic  temperature  profile  is  a 
satisfactory  approximation  for  calculating  transient  conduction  in  the 
solid.  Only  a small  error  is  thereby  caused  in  flamespreading  rate. 

(3)  Calculated  flamespreading  rate  is  sensitive  to  the  assumed 
ignition  criterion. 

(4)  Experimental  data  on  ignition,  transient  regression  rates, 
and  flamespreading  rates  are  too  scarce  to  form  a basis  for  definitive 
conclusions  on  the  correct  methods. 
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